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ABSTRACT 


The primary goal of this thesis was the development of an 
experimental technique, and supporting software, for the ac- 
quisition and analysis of unsteady velocity data generated by 
ayeeescillating airfoil. This research was in support of a 
major investigation of the compressibility effects on dynamic 
Stall. 

The experimental procedure involved schlieren flow visu- 
alization for comparison of steady and unsteady flow fields, 
and for determination of parameters for further study. Laser 
Doppler velocimetry was employed for obtaining velocity data 
in the airfoil wake. For unsteady data, the airfoil was 
eserlulated in pitch about its quarter chord. 

The data analysis produced wake profile plots represent- 
ing the flow field disturbed by the airfoil. Results were 


obtained for steady and unsteady conditions. 
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lL. INT RGBDUCTION 


a Dyer e STALL 
The possible benefits and consequences inherent sig 
eviamre °stal } have stimulated experimental investigations 


Pence Go Signer icanc progress im understanding this complex 


phenomenon. Much of the earlier work in this area was car- 
ried out by McCroskey using a two-dimensional Sir roa! 
Seestom TOSertivacing in sitch about its quarter-chord. The 


Results Geethis eftfort succeeded in eclarifying) the features 
Srecyhanmic Stallwup to a Mach number of 0.3 CRef. 13. 

Several Si@nificant distinctions between static and 
dynamic stall should be observed. Static stall is character- 
ized by massive flow separation and loss of lift coinciding 
Wot pee SOMe wr Cc rite a l angle of incidence. The aerodynamic 
forces on the airfoil vary uniquely, or nearly so, with angle 
ieoat tack. Dynamic stall, occuring during rapid increases 
in incidence, results ina delay of the associtated loss of 
lift until well beyond the static stall angle. Dynamic stall 
also tends to be considerably more persistant. The flow 
field does not immediately adjust to reductions in incidence 
angle waemnath static stall. CRef. 23 

Dynamic stall is characterized by the shedding and pas- 
Sa clemmecmmamvoertex over Che wpper siumrface Of an airfoil. The 


Preaquewlonmeor this vortex bewins shortly arter the incidence 


has exceeded that of static stall. A thin layer of reversed 
flow in the boundary layer, evolving into a vortex near ere 


leading edge, represents the initial phase of a dynamic stall 


eye le; As this disturbance moves across the upper surface cf 
the ale? oe be large negative pitching moments, SLENLficanme 
drag rise, and a-omonotonic incraasesin lift resulee The 
aerodynamic forces, reaching a Maximum as the vortex ate 


preaches the trailine edge, decrease dramatically with the 
sheddime vol em = were. Reattachment of the flow will then 
typically occur at some angle less than that corresponding tv. 


static stall. CRef. 2] 


Bund THES Poee@aAis 

The research project generating this thesis represents a 
small part of amajor investigation aimed at developing a 
thorough understanding of the effects of compressibility and 
Pitem oScrLiationm  ommdynmamitc sta wimw The primary goal of the 
research associated with this thesis was the development of 
an experimental procedure, and supporting software, for ene 
acquisition and analysis of unsteady flows with a laser dop- 
pler velocimetry system. The procedure and software devel- 
oped was to provide an initial analysis of possible areas of 
interest and a base upon which to build as the sophistication 
of the experimental methods increased. A successful attain- 
ment of this wzoal required ™ebtainine velocity data in the 
wake of an oscillating airfoil and Wake Profile ploees pro- 


duced from this data. 
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ieee EXPERIMENT AIC ONSTDERATIONS 


A. EXP ERD@MENTAL METHODS 
Sac e eee S Binecroduection a 1964, laser Doppler vel- 


cceimetry has proven to be a highly reliable means for obtain- 


ing measurements in complex turbulent flows. The method is 
Sapaele Sf ePperforming instantaneous, non-intrusive Ve loOcueey 
measurements Undem a Wide range ofr conditions. Being rela- 


tively independent of fluid properties, accurate measurements 
eS Pesst olen Ee Out tie requirement Of Calibration. LRef. 3] 

The basic principle behind laser Doppler velocimetry 
is the measurement of the velocity of particles travelling 
with the Polswael; These particles, when illuminated with a 
focused laser beam, become sources of scattered light. The 
Vopples shwreceof che scattered light is then weed to deter- 
Mine the veteettyegat the particles and, depending upon how 
well the particles are following the flow, the fluid velocity 
as well. 

AS a minimum, Six components are necessary for a com- 


plete LDV system: laser Source, transmitting optics, recelv- 


ing oOoptiies, PHotodetector , Signal processor, and data 
Processor The PURecerom- or the transmitting and Receiving 
Pee rccmare frOcuUSInNe sand light collecting, respectively. The 


EpHavcoGaeteoreer serves co convert lMeht Sstenal to an electrical 


Signal which is then used by the signal processor to provide 
a measurement of the frequency. The capability of deteacemmc 
reversed flows requires an additional component for frequency 
Shi mc ine - 

Due to the wide range of disciplines involved, Oop- 
taac Ss electronics, light scattering, and signal processing, 


the impleméentatiom of this method Can beudtfhiculcw a ec 35 


Consequently, the expertise of the operator can have a pro- 
found impact on the results of LDV measurements. Lt... ae 
therefore imperative, to an even greater degree than many 
other methods, that considerable attention be given to ap- 


plying correct measurement technique. 
e. Serlieres Plow Sse i ZaAcsoR 
Schlieren flow visualization systems provide a quali- 
tative analysis of the density gradient in compressible flow 
fields. The fundamental concept of this method is the vari- 


ation of the “refractive inde teradLont son mea fliid. Wives 


Gens Fey. 

A basic schlieren system includes a ijent  sovrcer 
transmitting and receiving mMigreors, a knife edge, and a 
focuSing screen. The transmitting mirror serves! €o directs 


parallel light beam through the test section and onto the 


rPecelvin@g.umirror. An image of €he light Setrece is then 
formed at the focal point of 2hiswmit tone The knife edge, 
placed in the plane of the light source image, cuts off pame 


af the transmitted ligeht and controls tne conv ega- eo 


image focused on the screen. Changes in test section density, 


therefore, peoguce  Varilatiomspor light intensity and provide 


a visual representation of the flow field. (CRef. 4) 


B. INSTRUMENTATION AND EQUIPMENT 
1. Laguipment Setup 

The experimental procedures for this research were 
conducted at the NASA Ames Fluid Mechanics Laboratory. A 
25-cm by 13-cm indraft wind tunnel with a slotted wall test 
section was used. Figure 2.1 shows the test segftion with 
Plenum chambers above and below the slotted walls. The 
adaptive wall design of this test section was developed fOr 


use in previous experiments. Although the technique was not 
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Figure 2.1 Wind tunnel test section. 


Ve 


used, the slotted walls were helpful in reducing Vital 
interference. 

The model for the experiment was a NAGA TOO 125a1- foi 
Withee.  —CHoOr al Vie valrrcil. attached to the test section 
Walls at quarter chord, was connected with mechanical linkage 


fouan -electricenetsar, AS ShONM Inthe ietiee sei. This variable 


(D 


speed motor provided ameans for oscillating the model about 
a Pitch axis anaes, obDtarniag unsteady data at various 
frequencies. 

An optical encoder, also driven by the motor, pro- 
vided 3,600 counts per revolution of 14-bit data on the phase 
angle of the oscillating airfoil. The characteristics of the 
encoder were absolute in that a whole word output, with a 
unique codespattern, represented each of the 3,600 discrete 
posileirense This code was derived from independent tracks on 
the encoder disk corresponding to individual photodetectors. 

Test section Mach number and velocity were obtained 
by using a 24-port Scanivalve with a differential pressure 
transducer. Calculation of these values required meas- 
urements from four Gf the ports, “ewo for Maem number computae 
tion and two for calibravien of ehe = ranceaucer. The linear 
relationship between the measured pressure and the corre- 
sponding voltage produced by the transducer was used fOr 
determining the pressures. In saeceomplishineg Gais, one pare 
of the Scanivalve, vented to the atmosphere, measured ambient 
pressure and a second port, tapped ta a mechanical vacuum 


Pump, measured a calibration pressure. A conversion’ )aceer 
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was then established by dividing the difference of the known 
pressures, Semospheric vanewcalibracionmweeyvecme Cifference in 


voltages produced by the transducer. 


Location of the LDV system probe volume, or beam 
crossing POinRe, was controlled by a three axis traversing 
mechanism. A Summit Taskmaster numerical control computer, 


eee tienally developed forscoentrolling movements of industrial 
machine tools, Was used for this task. The numerical control 
eomputer (NCC) was capable of manual operation, programmed 
eraratign. and, with software adaptation, operation from the 
csomputerswe@mkstation. VicmNe = Paey raed control inputs to two 
sets of electric motors, one for the transmitting optics and 
thHemmOcher fOoretne reecelvine optics. The necessity of using 
Serreremn.e traversing motors for the transmittine and receiv- 
inceope les placed a great demand on the synchronization of 
this system. Due to thewhigh senSpervicy of the LDV optics 
train, the synchronized movement was often inadequate for 
maintaining signal strength. 

Phase encoder data and digitized velocity data fram 
the LDV signal processors were channeled to a 3D-LDV Computer 
Interface CC hee The CI, multiplexing the inputs to a single 
Gigutal daitamwecnannel Towtput, provided 16-bit parallel output 
data to the computer data bus. mre Cl was capable of accept- 
ing up to six channels of digital data and from one to three 
event pulses from the signal processors. The multiplexing 
eycle could be initiated by random event pulse inputs or 


Svyrecnmnomized inputs, the latter being necessary for unsteady 


es. 


data acquisition: When operating in the synchronized mode, 
two or three pulses, as selected, occuring within a set time 
Ccoincidence time) would be required to initiate the multi- 
plex cycle and feed data to the computer. 
e- Laser Doppler Velocimeter Subsystems 

A two component LDV system was used to obtain the 
flow field velocity dee. Figures 2.2 and 2.3 show the gen- 
eral arrangement of the system components for the transmit- 
ting and receiving optics sides, respectively. 

A 5 Watt Argon ion laser, Spectra Physics Model 164, 


was used for the system. The unit was operated in the 


Seah eg 





Figure 2.2 LDV transmitting optics. 
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Figure 2.3 LDV receiving optics. 


multiline mode and produced a vertically polarized beam com- 
posed of all wavelengths inherent to this type of laser. 
Figure 2.4 illustrates the sequence of steps required 
of the transmitting optics. A set of mirrors was first used 
to lower the multiline beam to the level of the optical 
components. After passage EhEOuURH Aa ~“COllimator , used £oO 
ensure collimation of the beam and’ subsequent crossing of the 
beams at their waists, the beams pass through an attenuator. 
The function of this component was to reduce the beam intens- 
ity without inhibiting transmission of any wavelengths, a 
possible consequence of simply reducing the power output of 


the laser. Prior to separation of the wavelengths by the 


dispersion prism, it was necessary to horizontally polarize 


le: 
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Figure 2.4 Two-component transmitting optics. 
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Ene oeameoy passing it theoeugh a half-wave plate, er poliari- 
ZAtl On wGOorMator., invonrGcmeeomoobtain maximum light transmis- 
Son. The two stronger wavelengths, 514.5 nm (green) and 488 
mm €blue), were then guided by turning mirrors to the remain- 
ime OQPpwebeal components. 

Before “reaching  thetr respective Brage Cells, the 
green beam was split vertically and the blue was displaced 
to the center of the component assembly, Vemcicallya ipolar— 
ized, antics > leit a hotediz ont ally . The relationship between the 
Spaitting and pola@rizing= planes was driven by signal effi- 
eveone y. High §ienal—-to=-noise ratio requires Split beams to 
have equal intensities. The beams will be of approximately 
equal intensity when the polarity of the incident beam is 
perpendicular to the plane of the exiting beams CRef. 22), 
Therefore, the polarity of the vertically split beam was 
maintained asmnieg 2Zontal anamernat of the Norizentally split 
beam was rotated to vertical. As later discussed, this dif- 
ference in polarity was also necessary for separation of the 
BoaemenwaveCleneths by the) palarization beamsplitter in the re- 
Senwv ine Optics. 

Bragg Cells are acousto-optic devices used to gen- 
erate a frequency shift in the beam and allow for resolution 
of the 180 degree directional ambiguity. One beam of each of 
the two wavelengths was directed through a Brage Cell pra- 
tie ae sodas O Mize shi f t.. The unshifted beams were passed 
through glass rods to compensate for the path length differ- 


ence experienced by the shifted beam. Due to the nature of 


1B, 


the Bragg Cell, the frequency-shifted beam was deflected from 


L1ts inveial albvenmen-. Therefore, beam alignment prisms, or 


Tv 
a 


steering wedges, within the Bragg Cell modules served 
bring the beams back to parallel. 
A microscope ohiective, placed in the test section at 


the probe volume, was used tos observe the beams as adjust- 


ments were made to ensure the beams focused oO a commen 


Po tnien If the two probe volumes are not at the same location 
within the test section, then the measurements will be made 
at different perats= The necessity for all four beams to 


overlap at the focal point was especially prevalent for the 
acquisition of unsteady data due to the need for a Have 
coincidence data rate. 

Receiving optics are designed to image, to the great- 
est extent possible, the light scattered at the probe volume 
onto the photodetectors. Since the scattered light intensity 
is the greatest in the forward direction, a forward scattenm 
system is most often preferred for obtaining the greatest 
Signal-to-noise ratio. The system used for this analysis was 
an off-axis, forward scatteér configuration. Off-axis col= 
lection angles usually allow good noise isolation by admit- 
ting only light scattered at the sample volume i and improve 
measurement resolution by reducing the probe volume length. 
CRef. 3] 

The first three lenses in the receiving optics were 
used to collect scattered light from the probe volume and 


collimate the beam (Figure 2.5). The beam was then directed 
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to a Porat eaclion beamecelatcer for separation of the two 
SOmerS. The two components were then passed through line 
filters, for blocking ambient light and preventing cross 
talk, and to lenses used for focusing the beams onto the pin- 
Hele wae enaaelisemnof Chem raspeccive photomultiplier tubes. 
ie Pnietentpersbo lier tubes converted the light energy into 
electronic signals which were then passed to signal proces- 
Sing *efq{ulipment . Preamplifiers were located in close proxim- 
ity to the photomultiplier tubes and used to increase the 
Sena Level weerore the introduction of additional noise into 
Eine Cables. 

Slenal processing was accomplished with Macrodyne 
seunter Precessors. Counter processors operate by employing 
a crystal controlled ascillator to time over a selected num- 
ber wrer Cycles within the burst siegenal® Mmewelock count 1c 
initiated when the filtered signal exceeds a set threshold. 
The Doppler frequency is then determined by the time required 
te Count a specified number of cycles. The Macrodyne proces~- 
[Cuore lwoemautenemtne Clock and count over GS or 16 cycles, as 
selected. 

The Macrodyne processors use a periodicity check 
coupled with a three-level validation scheme to reject noise. 
The periodicity check uses two counters to make period meas- 
urements over a different number of cycles, 5S and 8 in this 
Ereceaure:. If the ratio of the two measurements is not with- 
in a specified percentage of 5/8, the signal is rejected. 


Tnhree-level validation requires a signal to exceed a positive 
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Figure 2.5 Data collection and processing. 
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wiresnala level, pass through zero, and then exceed a nega- 
tive threshold level before the next zero crossing. If these 
Semeletons — are not met Sem anyecyele the counter will resiet 
and immediately process a new Doppler burst. This method 
rrevents rejection of potentially good signals having noisy 
eve lies at the beginnine GE the burst. 

3 eG | Omens Syis tom 

The Scenhlieren system includéd two concave mirrors, 
Soe imeew ith a Pe eeine huet Ocal length, a Stropesece, 1c ee 2 te 
S Quvee e., a vertical knife edge, a focusing lens, and a polar- 
Pict inl SupPoaremtanregbtainims photegraphns of the flow field. 
The strobocscorfe, a type 3015 Strobrite, has a maximum output 
intensity of 2900 million beam candles and, using an internal 
oscillator, can provide flashing rates in excess of 1,000 per 
slecond. EVaeacdamelne a Lenseecoverswith a vertical golit to the 
Strobocecope A Slit source was obtained and the resolution of 
Price imase: Sigsgmaerrcantly improved. 

To obtain the unsteady photographs a digital compar- 
ator was used to trigger the stroboscope and provide the 
Peewee co Che phiatograrphic imaging on the’ pola@roid film. The 
Someacacormeerecetvod 14-bit phaseminput from =the optical en= 
coder and provided the triggering pulse when the 
instantaneous phase angle of the airfoil coincided with the 


angle specified €o the comparator. 
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aes DATA ACQUISITIOGN AUTOMATION 

A Digital VAXstation II minicomputer with graphics work- 
station was used to automate the data acquisition Drocesea 
The signal path aenaliustrated@em fFleure 2.0 start une from 
the probe volume, through the signal processing phases, and 
terminating within computer memory. The main program used in 
accomplishing this task, "LDV_ACQ’*, was originally developed 
to read data from the analog to digital converter. This was 
later modified to obtain pressure data for calculation of the 
test section velocity and call supporting subroutines for the 
acquisition and analysis of the velocity data. 

The subroutines assisting with various aspects of the 
acquisition process are listed below and included in” @ie 
appendix of this report. 

(1) HWTRVT - allows control of the traversing mechaniem 
from the computer terminal by sending character data 
to the interface. 

(2) HWWA - combined with an internally called subrou- 
tine, DAMUX, retrieves data from multiplexer, inter- 
prets phase and velocity information, creates data 
files, and generates histogram arrays and scaling 
parameters. 

(3)  PIJPLOT = provides graphical display of computed 
velocity values correspondinegr to airfoil phase 
angles. 

(4) PPLOT - similar to PT_PLOT But provides the addi. 
ional capability of allowing the user to vary the 


plotted velocity range. 


(5) HIST - plots a histogram of the percentage of data 
points which fall within specified velocity limits. 


IIIT. DATA ACQUISITION 


Aw LES SECTION VELOCITY 
The main program for the acquisition process, “EOY ACG % 


was responsible for the calculation of test section velocity. 


In combination with scanivalve positioning subroutines, the 
program acquired pressure transducer data from 7 ports. As 
previously discussed, measurements from 2 of the ports were 
used for calibration. Pressure data from 3 ports was’ dis- 


carded, and the other two provided the static pressure at the 
entrance to the test section and the total pressure measured 
at the nozzle inlet. Assuming isentropic, perfect gas condi- 
Tc lons 5 this supplied the necessary parameters for the Mach 
number computation by employing the following relationship: 
M = (((P,/P) = 1)" @/(y=1)"" 
with a value of 1.4 used for jy. The test section velocity 
was then calculated by: 
V=M xa 


where a, the sonic velocity, was obtained as follows: 


a=v yRT 
The temperature required for this calculation was read from a 
wall thermometer in the test cell and entered when prompted 
by the program. 
The velocity values were transferred to data files and 


later used to normalize the LDV velocity data. 


2° 


B. ( S€HoLIlLEREN PHOreeneaan. 

A qualitative analysis of the steady and unsteady flow 
fields was conducted at two Mach numbers, Os ane OS, and 
angles of incidence varying from 0 to 20 degrees. In aco- 
complishing this, test section Mach number was set by adjust- 
ing the wind tunnel throat and verified with Scanivalve 
pressures and calculations from ‘LDV_ACQ’. Final adjustments 
to the Stoboscope intensity and knife edge were then made to 


ensure optimum contrast. 


C. WARE VECOCG lt res 

With the desired wind tunnel conditions established and 
the probe volume in position for the data point, "HWWA’, a 
subroutine called by the main program, was employed for raw 


data acquisition and conversion to velocity and phase  infor- 


Mac ron. An internally called subroutine, ‘DAMUX’, acquired 4 
or 6 channels of data, as selected, from the LDV Computer 
[Inter falce El D- The CI was capable of transferring either 2, 


4, 6, or 8 channels of data, with the first two of the 
channels reserved for time and status information generated 
by the CI. During the acquisition of steady state velocities, 
4 channels were selected to provide the 2 velocity components 
in additven to time and status data, which were discarded 
Since they were not relevant for the experiment. For unsteady 
acquis terenm, 6 channels of data were transferred to accom- 


modate the additional data channel required for phase inputs 
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As previously mentioned, the signal processors determined 
the Doppler frequency by using the time required to count a 
specified number of cyles, Smrorethis analysis. This. raw 
data was transferred to the computer in 16 bits, with a 7106 
bit mantissa, a 4 bit exponent, and 2 nonessential bits used 
fom, controlling data flow. Due to the format of the time 
words, it was necessary to invert and reverse the bits before 
extracting the exponent and mantissa values. The signal 
period was then calculated as follows: 

ea—Gi sel (DIDS. ..B0 x 2 

where D9D8...D0 represented the binary mantissa and EXP was 
the range selected with the processor, used to divide the 
number of clock counts in the Doppler burst. 

The signal frequency, the inverse of the period, was ob- 


tained from: 


f = (32 x 1000)/(2EXP x mantissa) 


The signal frequency was composed of the Doppler frequen- 
cy, the 40 MHz shift frequency, and the downmixing frequency. 


This allowed the Doppler frequency to be deduced as follows: 


f=f—-f +f 


4 6 ehiff mis 


The 8 cycles counted within the Doppler signal physically 
represent 8 fringes crossed by a particle travelling with the 
flow. Consequently, the velocity component perpendicular to 
the fringe pattern may be calculated by the product of the 
Doppler frequency and the fringe spacing. The equation for 


the latter is: 
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6 = / [2 * sin(@/2)] 
where A is the wavelength of the beam and @ is the angle be- 
tween the intersecting beams. 
After obtaining a sufficient number of data points at a 


probe position (2,000 for unsteady flows), the velocity and 


phase data, as well as the probe position, test section 
velocity, and mixing frequencies, were stored in files for 
later analysis. The Taskmaster traverse controller was then 


used to move the probe volume and data collection SSR SEE <!. 
The data acquisition software included two routines’ for 
on-line performance checks of the data. “HIST*’ was used to 
provide a histogram of the percentage of data points which 
occurred within specified velocity limits. Figure 3.1 shows 
a sample plot provided by the subroutine. This data valida- 
tion method is most useful for steady flow conditions, in 
which case a high percentage of data concentrated around the 
estimated velocity would be a positive indicator. For Unl= 
steady analysis °PT_PLOT’ provided a graph of computed ve- 
locity values corresponding to airfoil phase angles. With 
this information available, the data quality could be 


easily assessed and a determination made regarding the ac- 


ceptance of these values. 
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Figure 3.1 Sample histogram. 
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IV. RESULTS AND DISCUSSION 


A> SCHLCIEREN PROTGGRAeHS 

The Schlieren photographs obtained from the analysis of 
the steady flow field were intended to provide information on 
the parameters corresponding to early flow separation and a 
reference for comparison with unsteady flows at, the same 
eonditions. Steady state separation was first observed at a 
Mach number of 0.5 and a 10 degree angle of incidence (Fig- 


ure 4.1). Based on this, the parameters for the unsteady 
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Figure 4.1 Steady state Schlieren photograph; 
mean a =10 degrees. 
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analysis were selected to be a Mach number of 0.5 and a phase 
angle corresponding to a 10 degree angle of incidence. Un- 
steady flow photographs were then obtained at oscillation 
frequencies of 5 and 25 Hz for comparison with each other and 
with wake profiles generated from LDV data (Figures 4.2 and 
Ghee 

From a comparison of the steady and unsteady photographs, 
it would appear that oscillating the airfoil reduces the size 


of the separated wake and, consequently, mitigates the effect 


of stalling. This is more apparent at the higher frequency 
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Figure 4.2 Unsteady flow Schlieren photograph; 
f=5 Hz, a =10 degrees. 
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Figure 4.3 Unsteady flow Schlieren photograph; 
f=25 Hz, a =10 degrees. 
oscillation, a logical result due to the greater effect of 


time unsteadiness at higher reduced frequencies. 


B. DATA ANALYSIS SOFTWARE 

The main program, "ORGIZE’, and its supporting subrou- 
tines were designed to provide wake profile plots from the 
data files generated during acquisition. To accomplish this, 
the data file corresponding to each point in the profile is 
input to the software. A series of routines first reorgan- 
izes the data files and discards statistically bad data 


points, those outside of two standard deviations from the 
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mean. Due to the oscillatory nature of the forcing function 
producing the unsteady data, a sinuSoidal curve fit is then 
introduced to allow evaluation of the first harmonic of the 
velocity distribution over the phase. APunique Sinusoid is 
generated fer each data point within the» profile. 

A minor degree of uncertainty is introduced into the ve- 
teeter cs bY thesg@ata interpolation and curve fitting rou- 
tines discussed above. However, unlike many other velocity 
measurement techniques, LDV data is obtained randomly due to 
the random arrival of the seed particles at the prove volume. 
Meet hout thisweurve fit, it would be difficult to egenerate 
velocity information at the same desired phase angle for each 
Gata WeOoint in the profile. 

During the generation of the wake profile the user is 


provided a comparison plot of mean values and the) sinusoid 


curve fit. The mean values are provided by ‘°BDPTS’ and com- 
puted within each 20 degree phase window, as illustrated in 
Figure 4.4. This was helpful in ascertaining the quality of 


the curve fit and, consequently, the quality of the wake pro- 
file results. 
The data analysis routines are summarized below and in= 
cluded in the appendix of this report. 
(1) ORGIZE - main program, calls supporting subroutines; 
sequentially organizes data to assist with subse- 


quent analysis. 


(2) PT_PLOT - provides graph of computed velocity values 
corresponding to phase. 


(3) BDPTS - computes mean and standard deviation of ve- 
locities occurring within 20 degree phase windows; 
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Figure 4.4 Comparison of mean and curve fit. 
discards data points outside of two standard devia- 
tions of the mean within each of these windows; pro- 
vides standard deviation plots of data. 

(4) AVG - averages data values occurring at same dis- 
crete phase. 

CS» INTERP - provides a linear interpolation between 
data, resulting in 180 equally spaced points; neces- 
sary for input to FOURIER. 

(6) FOURIER - computes coefficients for ae sinusoidal 


curve fit of the data.’adapted from June 1977 OCS 
Airfoil program). 


CZ» CVFIT - fits sinusoidal curve to velocity data Cu- 
tilizes FOURIER). 


a4 


CS) PLOT = provides wraphical comparison of curve folate 
with mean values generated by BDPTS. 


(9) PROPLOT - provides wake profile plots. 


Tae. EKOF TILES 


LDV data was obtained for steady and unsteady flows ata 


Mach number of 0.5 and a 10 degree angle of incidence. The 
steady profiles were collected at distances of 0.5, iO. eana 
23:0 inches behind the trailing edge of the airfoil. These 


values were non-dimensionalized by the chord length of the 
airfoil and are indicated as such on the profile plots. Due 
to time limitations, only one data profile was collected for 
the unsteady case, at two inches behind the trailing edge. 

The wake profile plots are presented in the following 
Pages of this section. A lack of data points in the lower 
part of the flow field will be apparent upon inspection of 
these results. This was due to a mechanical restriction on 
the traversing mechanism. 

The curve fit of the plotted values of the steady data 
WaS approximated and hand drawn to aid the reader in ob- 
serving the effects of the airfoil in the flow field. The 
software provides a spline interpolation fit between the 
data, as used in the unsteady plots. However, due to the 
amount of data scatter involved, this method was not used for 


these results. 


Sa) 


PROBE PUSTTION CIN. FRUM CENTER) 


HORIZONTAL VELOCITY PROFILE 





0.00 0.20 0.40 0.60 0.80 1.00 1220 


VELOCITY RATIO (U/UCINFIN) ) 


Figure 4.5 Steady wake profile. 
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Figure 4¢.6 Steady wake profile. 
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Figure 4.7 Steady wake profile. 
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Figure 4.8 Steady wake profile. 
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Figure 4.9 Steady wake profile. 
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Figure 4.10 Steady wake profile. 
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Figure 4.11 Unsteady wake profile. 
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Figure 4.12 Unsteady wake profile. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The purpose of the results obtained during the experi- 
mental phase of this thesis was to validate the technique and 
software developed. This was accomplished when unsteady 
data, generated by an oscillating airfoil, was successfully 
acquired, analyzed, and reduced to wake profile plots. 

Throughout the course of this investigation hardware 
problems have a posed a signigicant challenge to acquiring 
data within the existing time constraints. Consequently, the 
data obtained was accepted and included in this report with- 
out the opportunity to more carefully investigate regions of 
Significant data scatter and improve, or validate, the data 
plots. 

The results produced by these procedures confirm a sound 
technique, but also reveal several areas in need of addi- 
tional emphasis. At present, an interpolation subroutine is 
required to produce data at equally spaced intervals before 
generation of the curve fitting coefficients by *FOURIER®’. 
Combined with the inability to utilize a large number of data 
Pointe: 6 =< esoUum. this limitation will inherently introduce 
additional uncertainty into the results. A second software 
modification to be considered is the generation of the veloc- 
ity data at various phase angles simultaneously. Currently, 


the data must be processed for each separate phase angle 
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desired. This would eliminate the need for storage of large 
amounts of raw data and significantly reduce the processing 
time. By increasing the range of vertical movement of the 
traverse, the mechanical restriction that prevented the ac- 
quisition of data in the lower regions of the wake can be 
eliminated. 

Continuation of this study will be accomplished with a 
newly developed test section currently being installed. This 
test section will enable schlieren flow visualization and LDV 
data acquisition at locations over the airfoil which were 
previously blocked by the mechanical linkage used for sup- 
porting and oscillating the airfoil. This, combined with the 
modifications suggested above, will enable an in-depth analy- 
sis of wake and boundary layer characteristics associated 


with the dynamic stall phenomenon. 
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APPENDIX 


COMPUTER PROGRAMS 


PROGRAM LDV ACQ 


RARRAARARAHRARRARARARAANRARARAARARARARRAAAARARAARARARAAAARAARAARARHRAHRARARARHRARARARARAREES 


Adapted from LABSTAR Example3 1 by S. D. Hedrick in May, 1987 for reading 


the acquisition of unsteady velocity components with an LDV. 


ft 

t 

A-D data. Modified in August, 1987 to be used as the central program in * 
® 

f 

f 


RARARAKAKARAARARARARAREARARKRARARARARHRARARARARARAAARARARARARARAHRARARARAARARARARARAKRAHRARKRAAS 


INCLUDE ‘sysSlibrary:LIOSET. FOR’ 


C declare local variables 


INTEGER sys stat istatus ret by LIO calls 
INTEGER dev_id, clk_id !LIO device ID 
INTEGER data length tnumber of bytes of data read 


INTEGER i, npt, danpt 

INTEGER chn, nsteps, ier 

REAL rate, volt(0:23), press(0:23), corfac 

REAL sum, mean, calp, conv, amb, temp, mach, velo 
CHARACTER*1 ans 

CHARACTER*4 clock 

LOGICAL home 


C declare data buffers 


AAN 


AAADN 


10 


20 


INTEGER*2 raw _data(501) 
REAL VOLTAGES(501) 


OPEN(UNIT#13,FILE=’scani.dat’,STATUS=’NEW’ ) 


write(*,*) ‘ENTER A-D CHANNEL TO BE USED (0,1,8,9)’ 
read(*,*) chn 
werite(*,*) ‘ENTER CLOCK RATE (.01 ~ 4500 Hz)~-ALLOW FOR A’, 
‘ MAXIMUM OF 2 SECONDS TO ACQUIRE DESIRED SAMPLES’ 
read(*,*) rate 
if (rate .lt. .0l1 .or. rate .gt. 4500) then 
write(*,*) ‘RATE IS NOT WITHIN SPECIFIED LIMITS’ 
go to 10 
end if 
write(*,*)’ENTER NUMBER OF DATA POINTS TO BE COLLECTED’, 
‘ (500 max)’ - 


read(*,*) npt 

if (npt .lt. l°.or. net .gt.9 500) chen 
write(*,*) ‘NUMBER OR POINTS IS NOT WITHIN LIMITS’ 
go to 20 

end if 

npt = npt + 1 

dnpt = npt * 2 

write(*,*) ‘ENTER CLOCK DEVICE (KZA0O,KZBO) ’ 

read(*,’(A4)’) clock 


Attach the ADV11-D and set up for queued 1/0 


Sys stateLIOSATTACH(dev_id,’AZA0’, LIOSK QIO) 
IF(.NOT. sys Stat) CALL LIBSSIGNAL(tval(sys_ stat) ) 


Set up the ADV for synchronous transfer on selected channel 


Gain of 1 


Sys_stat = LIOSSET I (dev_id, LIOSK_SYNCH, 0) 
IF(.NOT. sys stat) CALL LIBSSIGNAL(tval(sys_ stat)) 
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ANMAANAN 


AAAN AAA 


AAA 


AaAN 


syS_ stat = LIOSSET_I (dev_id, LIOSK AD CHAN, 1, chn) 
IF(.NOT. sys stat) CALL LIBSSIGNAL(%val(sys_ stat) ) 


sys stat = LIOSSET I (dev_id, LIOSK AD GAIN, l, 1) 
IF(.NOT. sys_stat) CALL LIBSSIGNAL(%val(sys stat) ) 


Get a raw data buffer 
length of buffer and returned data length is in bytes 
Attach clock and set rate 


home = .true. 
nsteps = 0 
do 40 j-0,7 
call dalljllab(home,nsteps,ier) 
Sys _stat= LIOSATTACH(clk_id,clock,LIOSK QI0) 
IF(.NOT. sys_stat) CALL LIBSSIGNAL(%val(sys_ stat)) 
Sys _stat=LIOSSET Riclk_id,LIOS$K_CLK RATE, l,rate) 
IF(.NOT. sys stat) CALL LIBSSIGNAL(%val(sys stat)) 
Sys_stateLIOSSET I(cIk_id,LIOSK_TRIG, 1, LIOSK_ IMMEDIATE) 
IF(.NOT. sys stat) CALL LIBSSIGNAL(tval(sys_ stat)) 
sys stat = LIOSREAD(dev_id, raw_data,dnpt,data length, ) 
Irf(.NOT. sys_stat) CALL LIBSSIGNAL(%val(sys stat)) 
sys_ stat = LIOSSET I (dev_id, LIOSK TRIG, 2, LIOSK_CLK_ POINT, 
1 er id) 
IP(.NOT. sys stat) CALL LIBSSIGNAL(tval(sys_ stat)) 


Convert the raw data to voltages 


call LSPSFORMAT TRANSLATE ADC(raw data,voltages,npt, , , ) 


calculate the mean value of the voltages and detach the clock 


sum = 0.0 
do. 30" ie2,npe 
sum = sum + voltages(i) 

30 continue 

mean = sum / float(npt-1l) 

volt(j) = mean 

sys stat = LIOSDETACH(clk id, ) 

IF(.NOT. sys stat) CALL LIBSSIGNAL(tval(sys stat) ) 

40 continue 


detach from the A/D 


Ssys_stat = LIOSDETACH(dev_ id, ) 
IF(.NOT. sys stat) CALL LIBSSIGNAL(%val(sys_ stat) ) 


convert scanivalve voltages into units of pressure 


write(*,*) ‘ENTER CALIBRATION PRESSURE (in. of Hg)’ 

read(*,*) calp 

write(*,*) ‘ENTER REQUIRED CORRECTION FACTOR FOR CALIBRATION’, 
l ’ PRESSURE (in. of Hg)’ 

read(*,*) corfac 

calpecalp+corfac 

write(*,*) ‘ENTER ATMOSPHERIC PRESSURE (in. of Hg)’ 

read(*,*) amb 

press(0) = calp 
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C 


50 


55 


60 


press(1) © amb 
conv = (press(1l)-press(0)) / (volt(l)-volt(0)) 
do 50 je2,7 

press(j) # conv * volt(j) + calp - volt(0) * conv 
continue 
WRITE(*,*) VOLT(0),VOLT(1),VOLT(2),VOLT(7) 
write(l3755) 
format(2X,’PORT’,4X, 'PRESSURE’ ,//) 
do 60 i#0,7 

write(13,'(3X,12,5X,F9.4)’) i, press(i) 
continue 


C compute test section velocity 


C 


AY ANN 


70 


write(*,*)’ENTER AMBIENT TEMPERATURE (deg. F)’ 

read(*,*) temp 
machesgrt(abs(5.*((press(2)/press(7))**(.4/1.4)-1.))) 
veloemach*sgrt(1.4*53.3%32.174*( temp+459.67) ) 

write(*,*)’TEST SECTION MACH NUMBER AND VELOCITY-=--’ ,mach,velo 


write(*,*)’CONTINUE WITH LDV DATA ACQUISITION? (Y/N, DEFeyY)’ 
read(*,’(Al)’) ans 
if (ans .eq. ‘nm’ .or. ans .eq. ’N’) go to 70 


call subroutine to move LDV traversing mechanism and acquire data 


call HWTRVT(velo) 


stop 
end 
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AANANNANANANANANANNNANANAANANANNANNANDA 


ANNAANAANAANANANANANAANNANANANANNA 


SUBROUTINE HWWA(XPOS,ZPOS,VELO) 


ACRONYA: HardWare checkout - test of DRV11-WA 


PURPOSE: 
For LVIS, HWDVYB drives the acquisition of the LDV data. 


It is a test program for the YB driver. 
METHOD: 


BEGIN 

Install loadable drivers if necessary. 

Read the various instrument types in succession, skipping those which 
are not to be read. 

Remove loadable drivers if necessary. 

END 


PARAMETERS: 
ARG DIM LYPE 1/0/S2DESCRIPEION 


COMMONS USED: 
FILES USED: None 


ERROR HANDLING: 
IER »~ 1 means no “fatal” error. 


NOTES: 
LOCAL VARIABLES: 
VAR DIM TYPE DESCRIPTION 
EXTERNAL REFERENCES: 
NAME DESCRIPTION AND SOURCE 
DVWA Calls DRV11-WA device driver to acquire A/D data. 


STANDARDS VIOLATIONS: None. 
ENVIRONMENT: DEC VAX/VMS AND FORTRAN 77 


DEVELOPMENT HISTORY: 
DATE INITIALS DESCRIPTION 


07726743 TAL Adapted from DAPNT 
l2/7l2/7e> CLH Modified to print results in either 
octal or integer. 
02/25/87 GBG Modified for use with DRV11-WA and VMS 
05/06/87 CLH Modified to accept “new” DAMUX1. 
06/08/87 RRR Modified to subroutine + deleted CLH 5/87 Mod 
JUL-AUG,1987 SDH Modified to calculate velocities using variable 


mixing frequency as user input, interpret phase 
angle information generated by digital encoder, 
transfer test section velocity and probe posit- 
ion to data file, and generate histogram arrays 
and scaling parameters 


AUTHOR(S): Ted Lichtenstein, Informatics General Corp. 


jeRAARHKRAAHARAARARARKRAARARKRARARKRARARARARAARARARAARARHRERHKRERAARAEARARHKRARARARARARARARAREAL 


IMPLICIT NONE 
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Cc 


C 


ANNAN 


100 
1 
l 


200 


C 
220 


C 
C 


- 


PARAMETER NRAWMX=80000 
PARAMETER LUNTI #5 
PARAMETER LUNOUT#=6 
INTEGER* 2 IRAWDT(NRAWMX), IFORM, MUXTIME (16384) 
INTEGER*2 MUXSTS (16384), RAWDA( 20000), RAWDB( 20000) 
INTEGER* 2 RAWDC (20000) 
INTEGER STATUS, IER, PHAN 
INTEGER NCOLUM, IYESNO, NRAWIN,1,J 
INTEGER LIBSSUBX, NCHANS, NSAMPLES, INDEX, JNDEX 
INTEGER MRAWA, MRAWB ,MRAWC , MRAWD 
INTEGER IRAWDTA, IRAWDTB, IRAWDTC, IRAWDTD 
INTEGER IRAWDTE, IRAWDTF 
REAL FREQ(2),VEL(2),MIXA,MIXB,XPOS,2ZPOS,VELO 
LOGICAL DATSELCT (3) 
CHARACTER*10 INFILE 
CHARACTER®*1 ANS 
REAL * 4 YU(160),YV(120), YU1(160),YV1(120), YMAXU, YMAXV 
REAL* 4 SAVE 
INTEGER COUNT, XLOW, XHI 
REAL X(4000),y¥( 4000) 


INCLUDE ’(SSYSSRVNAM) ’ 


WRITE(*,*)’ ENTER DESIRED NAME OF DATA FILE--E.G., LDVXXX.DAT,’, 
’ WHERE THE "XXX" MAY BE USED AS A SERIAL ID FOR THE’, 
’ PROBE POSITION’ 

READ(*,'(A10)’) INFILE 

open (unitell, filesINFILE,status=’new’ ) 


CONTINUE 
Fill the raw data buffer with this bit pattern: 


"1010101010101010" 
This equals 125252(8), -21846(10). 


DO 200 I = 1, NRAWMX 
IRAWDT(1) = '125252'0 
CONTINUE 


DO I = 1, 16384 


MUXTIME (1) se) 12525200 
MUXSTS (1) = '125252'0 


END DO 


CONTINUE 
Prompt for no. of words to acquire: 


WRITE(*,*)’ ° 

WRITE(*,*)’Starting LDV acquisition’ 

WRITE(*,*)’ ’ 

WRITE(*,*)’ ENTER NUMBER OF CHANNELS TO BE USED’, 
’ (assumes non-oscillating airfoil if’, 
’ less than six)’ 

READ(*,*) NCHANS 


50 


224 WRITE A, 222) 
220 FORMAT (’ Enter number of samples: ’) 


READ (*, *) NSAMPLES 


IF (NSAMPLES .GT. 20000) THEN 
WRITE (*, 226) 


226 FORMAT (’ The number of samples must be 20,000 or less.’) 


GO TO 224 
END IF 


c 
230 CONTINUE 


AN 


aAaNN AANQANA 


aaAN 


* Prompt for output format of data: 


IFORM=#2 
UR () TFORM-LT-.1 .OR. IFORM.GT.2 ) GO TO 230 


Select the data to be returned. 


DATSELCT (1) @ .false. 
DATSELCT (2) @ .false. 
DATSELCT (3) @ .TRUE. 


Prompt for the mixing frequencies selected by user 


WRITE(*,*)’ ENTER MIXING FREQUENCIES FOR EACH COMPONENT - (blue, 
lgreen) - in MH2’ 
READ(*,*) MIXA,MIXB 


Read and display the data one sample at a time. 


IF (NCHANS .LT. 6) THEN 
WRITE(*,*)’ ENTER AIRFOIL PHASE ANGLE (tenths of degrees, 14)’ 
READ(*,*) PHAN 
DO 9 Iel,NSAMPLES 
CALL DAMUX(NCHANS,1,DATSELCT,MUXTIME,MUXSTS, IRAWDT, IER) 
RAWDA(I)@#IRAWDT( 1) 
RAWDB(I)#IRAWDT( 2) 
RAWDC (1) =PHAN 
9 CONTINUE 
GO TO 15 
END IF 
DO 10 Ieal,NSAMPLES 


CALL DAMUX (NCHANS, 1, DATSELCT, MUXTIME, MUXSTS, 
¢ IRAWDT, IER) 
RAWDA(I)#IRAWDT(1) 
RAWDB(1I)@IRAWDT( 2) 
RAWDC( I) =I RAWDT( 3) 
10 CONTINUE 


DO 19 1a1,NSAMPLES 
MRAWA=0 
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MRAWB=0 
MRAWC=0 
MRAWD=0 
MRAWA#=IIBITS(RAWDC(I),12,2) 
MRAWBeIIBITS(RAWDC(I),8,4) 
MRAWCeIIBITS(RAWDC(I),4,4) 
MRAWD=IIBITS( RAWDC(I),0,4) 
RAWDC( I) @MRAWA*1000+MRAWB*100+MRAWC* 1 0+MRAWD 
19 CONTINUE 
US: DO 20 I#1,NSAMPLES 
RAWDA(I)#INOT( RAWDA(I) ) 
RAWDB(I)=INOT( RAWDB(I) ) 
IRAWDTA=0 
IRAWDTB=0 
IRAWDTC#=0 
IRAWDTD=0 
IRAWDTE=0 
IRAWDTF#0 
IRAWDTA@IIBITS(RAWDA(I),10,4) 
IRAWDTBe-IIBITS( RAWDA(I),0,10) 
IRAWDTC#2**] RAWDTA* IRAWDTB 
IF (IRAWDTC .EQ. 0) THEN 
WRITE(*,*)’ IRAWDTC#0’ 
GO TO 20 
END IF 
FREQ(1)2@((32.*10.%*3.)/(PLOAT( IRAWDTC) ) )-40.+MIXA 
IRAWDTDeIIBITS( RAWDB(1),10,4) 
IRAWDTE@IIBITS( RAWDB(1I),0,10) 
- TRAWDTFe2**I RAWDTD* [IRAWDTE 
IF (IRAWDTF .£Q. 0) THEN 
WRITE(*,*)’ IRAWDTF#O’ 
f GO TO 20 
END IF 
FREQ(2)=((32.*10.**3.)/(FLOAT(IRAWDTF) ))-40.+MIXB 
vel(l)efreqg(2)*0.5145/(2.*sin(atan(13./(2.*%482.6)))) 
vel(2)#freq(1)*0.4880/(2.*sin(atan(13./(2.*%482.6)))) 
vel(l)evel(1)/.3048 
vel(2)evel(2)/.3048 


WRITE(11,700) vel(2),vel(1),RAWDC(I) 
20 CONTINUE 
WRITE(11,900) XPOS,Z2POS 
WRITE(11,950) VELO,NSAMPLES 
WRITE(11,975) MIXA,MIXB 


WRITE(*,*)’DO YOU DESIRE A HISTOGRAM OF THE VELOCITIES OBTAINED?’, 

1 ’ (Y/N - DEF@N)’ 

READ(*,’(Al)’) ANS 

IF (ANS .NE. ’y’ .AND. ANS .NE. ‘Y’) THEN 
GO TO 35 

END IF 

CLOSE( UNIT#11) 

OPEN(UNIT#1l] , PILE#INFILE,STATUS@='OLD’) 

DO 40 I#]1,NSAMPLES 
READ(11,800) X(I),Y¥(I) 

40 CONTINUE 

XLOW=-200 

XHI=-195 

DO 60 J-1,160 


Va 


AaAN 


ANAND 


COUNT=0 
DO 50 I=1,NSAMPLES 
IF (X(I) .GE. XLOW .AND. X(I) .LT. XHI) THEN 
COUNT#COUNT+1 
END IF 
50 CONTINUE 
YU(J)@( FLOAT( COUNT) /NSAMPLES)*100. 
XLOW=XLOW+5S 
XHI@XHI+¢5 
60 CONTINUE 


Calculating a scaling parameter for the horizontal velocity histogran 


DO 62 11,160 
YU1(I)eYU(I) 
62 CONTINUE 
DO 63 1e1,160 
DO 65 Jel,159 
IF (YU1L(J41) .LT. YUL(J)) THEN 
SAVE*#YU1(J) 
YU1(J)#YU1(J+1) 
YU1(J+1) SAVE 
END IF 
65 CONTINUE 
63 CONTINUE 
YMAXUeYU1(160)+10.0 


XLOWe- 300 
XHI=-295 - 
DO 80 Jel,120 
COUNT®20 
DO 70 Iel,NSAMPLES 
IF (Y(I) .GE. XLOW .AND. Y(I) .LT. XHI) THEN 
COUNT#COUNT+1 
END IF 
70 CONTINUE 
YV(J)@( FLOAT( COUNT) /NSAMPLES)*100. 
XLOW=XLOW+5 
XHI=XHI+5 
80 CONTINUE 


Calculating a scaling parameter for the vertical velocity histogram 


DO 82 Ie1,120 
YV1(T)eYV(TI) 
82 CONTINUE 
DO 85 J=1,119 
IF (YV1(Je1) .LT. YV1(J)) THEN 
SAVE#YV1(J) 
YV1(J)#Yv1(J+1) 
YV1(J+1) SAVE 
END IF 
85 CONTINUE 
83 CONTINUE 
YMAXV#YV1(120)+10.0 


CALL HIST( YU, YV, YMAXU, YMAXV) 


35 CLOSE(UNIT#11) 


nS 


WRITE(*,*)’DO YOU DESIRE A PLOT OF THE COMPUTED VELOCITY’, 


1 ' VALUES? (Y/N - DEFe#y)’ 
READ(*,’°(Al)’) ANS 
If (ANS .£EQ. ‘'N’ .OR. ANS .EQ. ‘n’) THEN 
GO TO 30 
END IF 


CALL PT PLOT(NSAMPLES, INFILE) 
WRITE(*,*)’DO YOU DESIRE TO VARY THE PLOTTED VELOCITY RANGE?’, 
l ' (Y/N--DEFe#yY)’ 
READ(*,’(A1l)’) ANS 
IF (ANS .EQ. ‘N’ .OR. ANS .EQ. ‘n’) GO TO 30 
CALL PPLOT(NSAMPLES, INFILE) 
* Again? 
30 WRITE(*,*)’ ° 
WRITE (LUNTI,*) ‘Again? (QO#NO, 1#YES) 
READ (LUNTI,*) IYESNO 
IF ( IYESNO.EQ.1 ) GO TO 100 


c 


RETURN 
* Focmats: 


700 FORMAT(5X,F9.4,1X,F9.4,8X,14) 

800 FORMAT(5X,F9.4,1X,F9. 4) 

900 FORMAT(//,2X,'XPOS#’ ,F6.2,2X,'ZPOS#=’ , F6.2) 

950 FORMAT(2X,’TEST SECTION VELOCITY#=’,F5.1,2X,15,’ SAMPLES’ ) 

975 FORMAT(2X,’BLUE AND GREEN MIXING FREQUENCIES--’,2F5.1,’(MHz)’) 


END 
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SUBROUTINE PT _PLOT(NSAMPLES, INFILE) 


a@eeee eae eaehakeananhanrakake ene aeaeneeheanreet eee e kane eehaaekhaneeaeeeeeateea ee ee kReae ete eaeeaenee 


* Written by S. D. Hedrick in July, 1987 to plot computed velocity values - 
* corresponding to airfoil phase angles. This subroutine is intended to « 
* provide the user with graphical information on acquired LDV velocities - 
* and assist in making a determination regarding the adequacy of the data - 
* obtained at the current probe position. * 
® « 
e e 


ARERR ARARARARRARERHARHRAERHRAKRAAEKRHRARHRARARARRHRERARARHRARKRARHARARHRARRARRARAERRERARARKRARAKRHRARARHR AREA 


IMPLICIT NONE 

INTEGER NSAMPLES,I 

REAL * 4 XVAL( 4000), YVALU( 4000) , YVALV( 4000) ,XCONTR( 4), YCONTR( 4) 

CHARACTER*10 INFILE 

CHARACTER*11 XLABEL 

CHARACTER*15 YLABEL 

CHARACTER*19 TITLE 

CHARACTER*1] DUMMY 

OPEN (UNIT#11,FILE#INFILE,STATUS=’OLD’ ) 

DO 10 Ie1,NSAMPLES 

READ(11,100) YVALU(1I),YVALV(I),XVAL(I) 
XVAL(I)@XVAL(1I)/10.0 
10 CONTINUE 
100 FORMAT(5X,F9.4,1X,F9.4,8X,F4.0) 

c 
C Set up the coordinate axis 
C 

XLABEL=’ PHASE ANGLE’ 

YLABEL=’VELOCITY (ft/s)’ : 

TITLE=’COMPUTED VELOCITIES’ 

XCONTR(1)=6.0 

XCONTR(2)=#0.0 y 

XCONTR(3)=360.0 

XCONTR(4)=#45.0 

YCONTR(1)=#6.0 

YCONTR( 2) =-200.0 

YCONTR(3)=600.0 

YCONTR(4)2#100.0 

CALL LGPSPLOT(1,’EXSY’,XVAL,YVALU,0,XLABEL,YLABEL, , , ,XCONTR, 

1 YCONTR, ,TITLE) 
Cc 
C Plot the individual data points on the axis 


CALL LGPSPOINT(1,XVAL,YVALU,.03,3,NSAMPLES, , , ) 
CALL LGPSPOINT(1,XVAL,YVALV,.03,4,NSAMPLES, , , ) 
CALL LGPSPUT_ TEXT(1,4.5,5.75,’* U-COMPONENT’, , ) 
CALL LGPSPUT_TEXT(1,4.5,5.5,'0 V-COMPONENT’, , ) 


WRITE(*,*) ‘Press <CR> to terminate plot’ 
READ(*,200) DUMMY 

200 FORMAT(A1) 
CALL LGPSTERMINATE PLOT(1,1) 
CLOSE(UNIT#=11) ia 
RETURN 
END 
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SUBROUTINE PPLOT(NSAMPLES, INFILE) 

IMPLICIT NONE 

INTEGER NSAMPLES, I 

REAL VMIN, VMAX, INTV 

REAL*4 xXVAL(4000),YVALU( 4000), YVALV(4000),xXCONTR(4),YCONTR(4) 

CHARACTER*10 INFILE 

CHARACTER*11 XLABEL 

CHARACTER*15 YLABEL 

CHARACTER*19 TITLE 

CHARACTER*1 DUMMY 

OPEN (UNITel1,FILE*INFILE,STATUS='OLD’ ) 

DO 10 I=1,NSAMPLES 
READ(11,100) YVALU(I),YVALV(I),XVAL(I) 
XVAL(I)=XVAL(I)/10.0 

CONTINUE 

FORMAT(S5X,F9.4,1X,P9.4,8X,F4.0) 


Set up the coordinate axis 


WRITE(*,*)’ENTER VELOCITY RANGE DESIRED (min,max)’ 

READ(*,*) VMIN, VMAX 

INTV=(VMAX-VMIN)/S. 

XLABEL@=’ PHASE ANGLE’ 

YLABEL=’VELOCITY (ft/s)’ 

TITLE=’COMPUTED VELOCITIES’ 

XCONTR(1)=#6.0 

XCONTR(2)<0.0 

XCONTR( 3)=360.0 

XCONTR(4)=45.0 

YCONTR(1)#6.0 

YCONTR(2)eVMIN 

YCONTR( 3) @VMAX 

YCONTR(4)©INTV 

CALL LGPSPLOT(1,’EXSY’ ,XVAL,YVALU,0,XLABEL,YLABEL, , , ,XCONTR, 
YCONTR, ,TITLE) 


Plot the individual data points on the axis 


200 


CALL LGPSPOINT(1,XVAL,YVALU,.03,3,NSAMPLES, , , ) 

CALL LGPSPOINT(1,XVAL,YVALV,.03,4,NSAMPLES, , , ) 

CALL LGPSPUT TEXT(1,.05,.3,’ASTERISK=U COMPONENT’, , ) 
CALL LGPSPUT TEXT(1,.05,.1,’CIRCLE*V COMPONENT’, , ) 


WRITE(*,*) ‘Press <CR> to terminate plot’ 
READ(*,200) DUMMY 

FORMAT(A1 ) 

CALL LGPSTERMINATE PLOT(1,1) 
CLOSE(UNIT#11) 

RETURN 

END 
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SUBROUTINE HIST( YU, YV,YMAXU, YMAXV) 


RAARARAARARAARRARARAAARARARARARARARAEAARHERAARARAAARAARARRRAR ARAN 


Written by S. D. Hedrick in August, 1987 to plot a histogram of the per- 
centage of data points which fall within specified velocity limits. The 


ing parameters for the vertical axis and data arrays are computed within 


a preceding subroutine (HWWA). 


2 

® 

® 

* subroutine produces histograms for each of two velocity components. 
® 

® 

® 

® 


Raeaeeazeekaaeeanee eee aeeee ate e eae ee hee eee ee enkhaeee eee eee eaeeeeeenaktee ee eee tea eee nee eee 


C 


IMPLICIT NONE 

REAL*4 YU(160),YV(120),XCONTRU(4),XCONTRV(4),YCONTRU( 4) 
REAL*4 YMAXU, YMAXV, YCONTRV( 4) 

REAL XLOW,XHI,XLOWU(160),XHIU(160) ,XLOWV(120),XHIV(120) 
INTEGER I 

CHARACTER*29 UTITLE 

CHARACTER*27 VTITLE 

CHARACTER*10 XLABEL 

CHARACTER*18 YLABEL 

CHARACTER*1 DUMMY 


C Setting up the coordinate axis 


c 


C 


XLABEL#’VELOCITIES’ 

YLABEL#=’ PERCENTAGE OF DATA’ 

UTITLE=’HORIZONTAL VELOCITY HISTOGRAM’ 

VTITLE='VERTICAL VELOCITY HISTOGRAM’ 

XCONTRU(1)#=#6.0 

XCONTRU( 2) #-200.0 2 

XCONTRU( 3)=#600.0 

XCONTRU(4)2#100.0 

XCONTRV(1)=#6.0 

XCONTRV(2)=#2-300.0 

XCONTRV( 3)=#300.0 

XCONTRV(4)2#100.0 

YCONTRU(1)=#6.0 

YCONTRU(2)2#0.0 

YCONTRU (3) =YMAXU 

YCONTRU(4)#10.0 

YCONTRV(1)#=#6.0 

YCONTRV(2)#0.0 

YCONTRV( 3) =#YMAXV 

YCONTRV(4)=#10.0 

CALL LGPSPLOT(1,’IXSY’, ,YU,0,XLABEL,YLABEL, , , ,XCONTRU, 
YCONTRU, ,UTITLE) 


C Setting up arrays for horizontal and vertical histogram bar parameters 


C 


10 


XLOW=-200.0 

XHI#-195.0 

DO 10 Ie1,160 
XLOWU( I) @XLOW 
XHIU(I)#XHI 
XLOWeXLOW+S.0 
XHI@XHI+5.0 

CONTINUE 

XLOW=-300.0 

XHIe-295.0 

DO 20 121,120 
XLOWV(1I)#XLOW 


S7 


Scal- 


& 
® 
* 
® 
® 
® 
® 
® 


nm NAAN 


Anan 


XHIV(I)*XHI 

XLOW=XLOW+5. 0 

XHI@XHI+5.0 
20 CONTINUE 


Plotting the horizontal velocity histogram 
CALL LGPSHIST(1,XLOWU,XHIU,YU,160, , ,2) 


WRITE(*,*) ‘PRESS <CR> TO TERMINATE PLOT’ 
READ(*,‘’(Al)’) DUMMY 
CALL LGPSTERMINATE_PLOT(1,1) 


Repeating the procedures to plot the vertical velocity histograas 


CALL LGPSPLOT(1,’IXSY’, ,YV,0,XLABEL,YLABEL, , , ,XCONTRV, 
1 YCONTRV, ,VTITLE) 
CALL LGPSHIST(1,XLOWV,XHIV,YV,120, , ,2) 


WRITE(*,*) ‘PRESS <CR> TO TERMINATE PLOT’ 
READ(*,’(Al)') DUMMY 

CALL LGPSTERMINATE_PLOT(1,1) 

RETURN 

END 
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PROGRAM ORGIZE 


eeaetaeetaeetkheeetae teeta eeeeeetaeaeeeteeteeeeeeeeeettheeeaeeeetaeeeeeateteheeheeeeeeeees 


e 
« 
2 
« 
& 
2 
« 
& 
® 
® 
® 


Written by S. D. Hedrick in July, 1987 to sequentially organize data files 
to obtain a least value of phase angle (third column of data), and associ- 
ated velocities, at top of file. This rearrangement of data was desired 
to assist with subsequent manipulation and analysis. In September, 1987 
this routine was incorporated as the central program in the processing and 


analysis of unsteady LDV velocity data generated by an oscillating airfoil. 
The program calls assisting subroutines to disgard bad data eames and pro- 


vide wake profile plots. 


Seek eecaeetekeeeteeeeee teehee ea hee teaetekeeatheae thee thee etehehataktteeatektaeke tee tehathathate 


IMPLICIT NONE 
REAL X(4000),Y(4000),XPOS,2ZPOS,VELO,UVL,VVL,UNDIM( 30) ,VNDIM( 30) 
REAL 2P0(30),SAV1,SAV2,MENU,MENV, PREJ,MEANU( 18) ,MEANV(18) 
REAL MPHAVG( 18) 
INTEGER 2(4000),NSAMPLES, PHASA,I,J,K,N,M,SAVE,REJ,STDY 
CHARACTER*10 INFILE 
CHARACTER*1 ANS 
STDY=#0 
WRITE(*,*)’IS ANALYSIS FOR NON-OSCILLATING AIRFOIL? (Y/N)’ 
READ(*,’(Al)’) ANS 
IF (ANS .EQ. ’Y’ .OR. ANS .EQ. ‘y’) STDYe#]1 
WRITE(*,*)’ ENTER NUMBER OF DATA ACQUISITION POINTS IN PROFILE’, 
| ‘ (30 maximum)’ 
READ(*,*) N 
WRITE(*,*)’ ENTER PHASE ANGLE DESIRED FOR ANALYSIS (degrees) ’ 
READ(*,*) PHASA 
DO 70 Kel,N 
WRITE(*,*)’ ENTER NAME OF INPUT DATA FILE (LDVXXX.DAT)’ 
READ(*,’(A10)’) INFILE 
WRITE(*,*) ’ENTER NUMBER OF SAMPLES IN DATA FILE’ 
READ(*,*) NSAMPLES 
OPEN(UNIT#11,FILE#=INFILE,STATUS#=’OLD’ ) 
DO 10 Ie#l1,NSAMPLES 
READ(11,100) X(1I),Y¥(I),2(T) 


10 CONTINUE 


READ(11,200) XPOS,ZPOS 
READ(11,300) VELO 
mrMesustOY 7EO. 1) THEN 
CLOSE(UNIT#11) 
CALL STEADY (NSAMPLES,INFILE,MENU, MENV) 
UNDIM(K)=«MENU/VELO 
VNDIM( K)=#=MENV/VELO 
GO TO 80 
END IF 
M=eNSAMPLES-1 
DO 20 J=1,NSAMPLES 
DO 30 I1#1,M 
BE (20isi) {LT. Z(1)) THEN 
SAVE=2(T) 
SAV1=«X(T) 
SAV2#Y(T) 
Z(I)e2(T+1) 
X( I) eX(T+1) 
Y(I)#e¥(T1) 
Z(3+1)°SAVE 
X( I+] )=SAV1 
Y¥(I¢1)=SAV2 


ee 


« 

2 

@ 
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100 
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END IF 

CONTINUE 
CONTINUE 
CLOSE(UNIT#11) 
OPEN(UNITe#15, FILE*INFILE,STATUS='’OLD’ ) 
DO 40 I=#1,NSAMPLES 

WRITE(15, 200) =X(1),Y(1),2(1) 
CONTINUE 
CLOSE(UNITe#15) 
WRITE(*,*)’DO YOU DESIRE A VELOCITY/PHASE PLOT OF 


n ’ DATA? (Y/N--DEFey)’ 


READ(*,’(Al)’) ANS 

IF (ANS .EQ. ‘'N’ .OR. ANS .EQ. ‘’n’) GO TO 50 

CALL PT_PLOT(NSAMPLES, INFILE) 

CALL BDPTS(NSAMPLES, INFILE,MENU,MENV,REJ, PREJ,MEA: 


1 MPHAVG) 


CALL AVG(NSAMPLES,INFILE) 

WRITE(*,*)’DO YOU DESIRE INTERPOLATION OF THE DAT: 
READ(*,’(Al)’) ANS 

IF (ANS .EQ. ‘'N’ .OR. ANS .EQ. ‘n’) GO TO 55 

CALL INTERP(NSAMPLES, INFILE) 

CALL CVFIT(NSAMPLES, INFILE, PHASA, UVL, VVL, MENU, MEN\ 


1 MEANU, MEANV, MPHAVG) 


UNDIM(K)*#UVL/VELO 
VNDIM(K)@#VVL/VELO 
ZPO(K)=#ZPOS 

CONTINUE 

CALL PROPLOT(UNDIM,VNDIM,2PO,XPOS,N) 


- 


Format statements 


FORMAT (5X,F9.4,1X,F9.4,8X,14) 
FORMAT(//,7X%,f6.2, /X,1D6. 2} 
FORMAT (24X,F5.1) 

STOP 

END 
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40 


SUBROUTINE BDPTS(NSAMPLES, INFILE,MENU,MENV,REJ,PREJ, 
1 MEANU, MEANV, MPHAVG ) 
REAL X(4000),Y(4000),MEANU(18),MEANV(18),STDVU(18) 
REAL STDVV(18),MPHAVG(18),MINU,MINV, MAXU, MAXV, MENU, MENV 
INTEGER 2(4000),PHASEA,PHASEB, SUM, COUNT,NSAMPLES,NSAMP,REJ 
REAL*4 XCONTR(4),YCONTR(4),PREJ 
CHARACTER*1 ANS,DUMMY,DUMMI 
CHARACTER*10 INFILE 
CHARACTER*11 XLABEL 
CHARACTER*21 YLABEU 
CHARACTER*]19 YLABEV 
CHARACTER*18 TITLE 
OPEN(UNIT#15,FILE#INFILE,STATUS#=’OLD’ ) 
NSAMP#NSAMPLES 
DO 10 I#l,NSAMPLES 
REND(1S,100) X(r),Y(I),2(T) 
CONTINUE 
FORMAT(5X,F9.4,1X%,F9.4,8X,14) 
PHASEA#0 
PHASEB#200 
DO 20 #1,18 
SUM=0 
SUMUs0. 
SUMV=Q. 
PHAVGe0. 
DO 30 J=al,NSAMPLES 
IF (Z(J) .LE. PHASEB .AND. 2(J) .GE. PHASEA) THEN 
SUMU=SUMU+X( J) 
SUMV=SUMV+Y(J) 
PHAVGe PHAVG+ ( FLOAT(2(J))/10.) 
SUM=SUM+1 
END IF 
CONTINUE 
IF (SUM .EQ. 0) THEN 
MEANU(I)#0. 
MEANV(I)#0. 
MPHAVG(I)#=(FLOAT( PHASEA)+100.)/10. 
VARUsO. 
VARVe(. 
GO TO 40 
END IF 
MEANU(I)=SUMU/FLOAT( SUM) 
MEANV(I)=SUMV/FLOAT( SUM) 
MPHAVG(I)=PHAVG/FLOAT( SUM) 
VARU=0. 
VARV=0. 
DO 40 J=l1,NSAMPLES 
IF (Z(J) .LE. PHASEB .AND. 2(J) .GE. PHASEA) THEN 
Al=X(J)-MEANU(I) 
A2=Y(J)-MEANV(I) 
BleAl**2/FLOAT( SUM) 
B2=A2**2/FLOAT( SUM) 
VARU*VARU+B1 
VARV=VARV+B2 
END IP 
CONTINUE 
STDVU(1I)*eSQRT(VARU) 
STDVV(I)*eSQRT( VARV) 
PHAS EA@ PHAS EA+200 
PHAS EB=PHASEB+200 


oe 
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CONTINUE 
WRITE(*,*)’DO YOU DESIRE A STANDARD DEVIATION PLOT? (Y/N - DEF#Y)’ 
READ(*,200) ANS 
FORMAT(A1 ) 
IF (ANS .£Q. 'N’ .OR. ANS .EQ. ‘n’) THEN 
GO TO 50 
END IP 


Setting up the coordinate axis for the plot 


XLABEL@#’ PHASE ANGLE’ 

YLABEU=}/HORIZONTAL VELOCITIES’ 

YLABEV#’VERTICAL VELOCITIES’ 

TITLE@#’ STANDARD DEVIATION’ 

XCONTR(1)#6.0 

XCONTR(2)20.0 

XCONTR(3)=360.0 

XCONTR(4)#45.0 

YCONTR(1)=#6.0 

YCONTR(2)=#-200.0 

YCONTR(3)=#600.0 

YCONTR(4)#100.0 

CALL LGPSPLOT(1,’EXSY’,MPHAVG,MEANU,0,XLABEL,YLABEU, , , ,XCONTR, 
YCONTR, ,TITLE) 


Plotting the horizontal velocity data points on the axis 


CALL LGPSPOINT(1,MPHAVG,MEANU, .03,4,18, , , ) 


Plotting the horizontal standard deviation 


45 


CALL LGPSSTNDEV(1,MPHAVG,MEANU,STDVU,18, ) : 


WRITE(*,*) ‘PRESS <CR> TO TERMINATE PLOT’ 
READ(*,200) DUMMY 
CALL LGPSTERMINATE_ PLOT(1,1) 


Plotting the vertical standard deviation 


47 


50 


CALL LGPS$SPLOT(1,’EXSY’ ,MPHAVG, MEANV,0,XLABEL,YLABEV, , , ,XCONTR, 
YCONTR, ,TITLE) 

CALL LGPSPOINT(1,MPHAVG,MEANV,.03,4,18, , , ) 

CALL LGPSSTNDEV(1,MPHAVG,MEANV,STDVV,18, ) 

WRITE(*,*) ’PRESS <CR> TO TERMINATE PLOT’ 

READ(*,200) DUMMI 

CALL LGPSTERMINATE PLOT(1,1) 


COUNT#1 
PHASEA#0 
PHASEB=200 
DO 60 1r#1,18 
MINU@eMEANU(I)-2.0*STDVU(I ) 
MINVeMEANV(I)-2.0*STDOVV(I) 
MAXUeMEANU(I)+2.0*STDVU(I) 
MAXVeMEANV(1I)+2.0*STDVV(I) 
DO 70 Jel,NSAMPLES 
IF (2(J) .LE. PHASEB .AND. 2(J) .GE. PHASEA) THEN 
IF (X(J) .LT. MINU .OR. X(J) .GT. MAXU) THEN 
COUNT#=COUNT+1 
IF (J .E£Q. 1) GO TO 70 
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80 


85 


70 


60 


a0 


95 


97 


98 
96 


DO 80 K=0,J-2 
X(J-K)@X(J-K-1) 
Y¥(J-K)eY(J-K-1) 
Z(J-K)e2(J-K-1) 

CONTINUE 

GO TO 70 

END IF 
IF (Y¥(J) .LT. MINV .OR. Y(J) .GT. MAXV) THEN 

COUNT=COUNT+1 

tet) 7.6Q. 1) GO TO 770 

DO 85 K#0,J-2 
X(J-K)=X(J-K-1) 
Y¥(J-K)#Y¥(J-K-1) 
Z(J-K)#Z(J-K-1) 

CONTINUE 

END IF 
END IF 
CONTINUE 
PHASEA=PHASEA+200 
PHASEB=PHASEB+200 
CONTINUE 
CLOSE(UNIT#15) 
OPEN( UNIT#12,FILE*=INFILE,STATUS#=‘OLD’ ) 
DO 90 I*#COUNT,NSAMPLES 
WRITE( 12,100) 4(1)7x(1),2(1) 
CONTINUE 
NSAMPLES=NSAMPLES+1-COUNT 
CLOSE(UNIT#12) 
OPEN(UNIT#13,FILE#INFILE,STATUS#=’OLD’ ) 
Suus=0. 
SUVeO. 
DO 95 I=#1,NSAMPLES 
READ( 13,200) X( 1), Y(1),2(1) 
SUUs#SUU+X( 1) 
SUV=SUV+Y(TI) 
CONTINUE 
MENU@#SUU/FLOAT(NSAMPLES ) 
MENV=SUV/FLOAT(NSAMPLES ) 
CLOSE(UNIT#13) 
REJ@eNSAMP-NSAMPLES 
PREJ=#=(1.-(FLOAT(NSAMPLES) /FLOAT(NSAMP)))*100. 
PHASEA#=0 
PHASEB#200 
DO 96 181,18 
SUM=0 
SUMU=0. 
SUMVe0. 
PHAVG=0. 
DO 97 Jal,NSAMPLES 
IF (Z(J) .LE. PHASEB .AND. Z(J) .GE. PHASEA) THEN 
SUMU*#SUMU+X(J) 
SUMV=SUMV+Y(J) 
PHAVG*PHAVG+( FLOAT(Z(J))/10.) 
SUM=SUM+1 
END IF 
CONTINUE 
IF (SUM .EQ. 0) THEN 
MEANU(1)=0. 
MEANV(I)=0. 
MPHAVG(1)#=( FLOAT(PHASEA)+100.)/10. 
GO TO 98 
END IF 
MEANU( 1) *SUMU/FLOAT( SUM) 
MEANV( 1) =SUMV/FLOAT( SUM) 
MPHAVG(1)=*PHAVG/FLOAT( SUM) 
PHAS EA=PHASEA+200 
PHASEB=PHASEB+200 
CONTINUE 
RETURN 
END 
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SUBROUTINE AVG(NSAMPLES, INFILE) 
REAL X(4000),Y(4000) 
INTEGER 2(4000),COUNT,SUM,NSAMPLES 
CHARACTER*10 INFILE 
COUNT#]1 
M=NSAMPLES-1 
OPEN(UNIT#12,FILE=INFILE,STATUS#=’OLD’ ) 
DO 10 I=#1,NSAMPLES 
READ(12,100) X(I),Y¥(1I),2(1T) 
CONTINUE 
FORMAT(5X,F9.4,1X,F9.4,8X,14) 
DO 20 I=#1,M 
IF (2(I41) .£EQ. 2(1) . AND. Z(I+2) .NE. Z(1)) THEN 
COUNT=COUNT+1 
X(I4+1)=e(X(1+1)4X(I))/2.0 
Y(I+#l)e( y¥(I4+1)4+Y(1I))/2.0 
DO 30 J#0,I-2 
X(I=-J)eX(I-J-1) 
Y¥(T-J )#Y(I-J-1) 
Z(I-J )eZ(I-J-1) 
CONTINUE 
END IF 
IF (Z(1+2) .£Q. Z2(1)) THEN 
SUM=1 
DO 25 Je2,NSAMPLES-1 
IF (Z(14J}) .EO. Z(1))9 THEN 
SUM=SUM+1 
END IF 
CONTINUE 
DO 27? J=1,SUM 
X(I+SUM) =X(1+SUM)+X(1I+SUM-J) 
Y¥(I1+SUM)#Y(1+SUM)+Y(1I+SUM-J) 
CONTINUE 
X(I+SUM) =X(1I+SUM) /(SUM+1 ) 
Y(I+SUM)#Y(1I+SUM)/(SUM+1) 
DO 28 Jel,I-1 
X(I+SUM-J)eX(I-J) 
Y(1+SUM-J)#Y¥(I-J) 
2(1+SUM-J )eZ(I-J) 
CONTINUE 
COUNT=COUNT+SUM 
END IF 
CONTINUE 
CLOSE( UNI T#12) 
OPEN(UNIT#14,FILE#INFILE,STATUS#=’OLD’ ) 
DO 40 IeCOUNT,NSAMPLES 
WRITE(14,100) X(1I),Y(I),2(T) 
CONTINUE 
NSAMPLES#=NSAMPLES+1-COUNT 
CLOSE(UNIT#14) 
RETURN 
END 
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SUBROUTINE INTERP(NSAMPLES, INFILE) 
IMPLICIT NONE 
CHARACTER*10 INFILE 
REAL X(3600),yY(3600) 
INTEGER 2(3600),NSAMPLES,1,J,COUNT,NSA, FAC 
CHARACTER*1 ANS 
OPEN(UNIT#11,FILE*INFILE,STATUS=’OLD’ ) 
DO 10 I=#1]1,NSAMPLES 

RenO( 11,100) X(I),Y(1),.zZ(1} 
CONTINUE 
FORMAT(5X,F9.4,1X,F9.4,8X,14) 
COUNT#=0 
NSA*NSAMPLES 
DO 20 11,3600 

IF (COUNT .GE. NSA) THEN 

Z(I)eI-1 


RCI )=((X(1)-X( 1-1) )/(2(1)-2(1-1))) * (201) -20 1-1) )+x(1I-1) 
Y(T)o( (YL) -Y¥( 1-1) ) 70201) -2( 1-1) ))*(2(1)-2( 1-1) )e¥( 1-1) 
NSAMPLES=NSAMPLES+1 
GO TO 20 

END IF 


IF (2(1) .£Q. I-1) COUNT=sCOUNT+1 
IF (Z(I) .NE. I-1) THEN 
NSAMPLES*NSAMPLES+1 
DO 30 J=0,NSAMPLES-I-1 
X(NSAMPLES-J)=X(NSAMPLES-J-1) 
Y(NSAMPLES-J)#Y(NSAMPLES-J-1) 
Z(NSAMPLES-J)}©Z%(NSAMPLES-J-1) 
- CONTINUE 
Z(I)=I-1 
Piet) BOwe ll) THEN 
X(T )e((X(1+1)-X(NSAMPLES) )/(Z2(1+1)-Z2(NSAMPLES) ))*(Z2(1)- 
Z(NSAMPLES ) )+X(NSAMPLES ) 
Y(T)|((Y(1I+1)-Y(NSAMPLES) )/(Z2(1+1)-Z(NSAMPLES)))*(Z(1I)- 
Z(NSAMPLES) )+Y(NSAMPLES) 


GO TO 20 
END IF 
gla er ere te MCT — 1) 7 (2(14+1)<-2(1=-1)))*(2(1)-2(1=1) }+x(i-1) 
Vier rey —y(I-)))/(2( 141 )-Z2(1=1)))*(2Z{1)-Z2(1=-1))+Y( 1-1) 
END IF 
CONTINUE 


CLOSE(UNIT#11) 

WRITE(*,*)’DO YOU DESIRE TO REDUCE THE NUMBER OF INTERPOLATED’, 
1 ’ DATA POINTS BY AN INTEGER FACTOR OF 3600? (Y/N)’ 
READ(*,’(A1l)’) ANS 
OPEN(UNIT#12,FILE*#INFILE,STATUS=’OLD’ ) 

IF (ANS .EQ. ‘Y’ .OR. ANS .EQ. ‘y’) THEN 
WRITE(*,*)’ ENTER INTEGER FACTOR’ 
READ(*,*) FAC 

FAC=20 
DO 34 IseFAC,NSAMPLES, FAC 
WRITE(12,100) X(1I),Y¥(I),2(1) 
CONTINUE 
NSAMPLES=3600/FAC 
GO TO 40 

END IF 

DO 40 Iel,NSAMPLES 
WRITE(12,100) X(I),¥(I),2(1) 

CONTINUE 

CLOSE(UNIT#12) 
RETURN 
END 
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SUBROUTINE FOURIER(F,M,A,B,N) 
DIMENSION F(0:M),FA(0:300),FB(0:300) 
DATA PI/3.14159265/ 


He2*PI/M 

TYPES eed 

J=0 

DO 2 I#0,M 

Ca4 

IF (1 .80. J)g¢=z 

IF (I .EQ. J) JaJe2 
FA(I)=F(1I)*COS(N*H*TI) 
FB(I)eF(1I)*SIN(N*H*T ) 

IF (I .EQ. M) T3eFA(T) 

IF (1 SEQ. M)RITISFS CI) 
IF (I .£Q. M-1) T2=FA(T) 
IF (I .EQ. M-1) TT2=FB(1) 
IF (I .EQ. M-2) TleFA(TI) 
IF (I .EQ. M-2) TT1=FB(I) 
IF (I .EQ. M-3) TO#FA(TI) 
IF (I .EQ. M-3) TTO#=FB(I) 
IF (I) EQ. 0) sGesto 2 


FA(I)*#FA(I-1)+C*FA(T) 
FB(I)@#FB(I-1)+C*FB(TI) 

CONTINUE 

AeH/3./PI*(FA(M)-T3) 
BeH/3./P1I*(FB(M)-TT3) 

IF (J .NE. M+#1) GO TO 3 
AwH/3./PI*(FA(M)-4.*T3-2.*T2-4.*T1-TO) 
AwA+3.*H/8./PI*(TO+3.*T14+3.*T2+T3) 
BeH/3./PI*(FB(M)-4.*TT3-2.*TT2-4.*TT1-TTO) 
BeB+3.*H/8./PI*(TT0+3.*TT1+3.*TT2+TT3) 
CONTINUE ‘ 
RETURN 

END 
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SUBROUTINE CVFIT(NSAMPLES,INFILE, PHASA, UVL, VVL,MENU,MENV,REJ,PREJ, 
1 MEANU, MEANV , MPHAVG ) 


ez#eenweeeettheocnee ea teeteaehteetee teehee eneekteeteaeneeneeteeeeeeeeeee eee eene eee eee ete eea ee eeee 


e 
* Writtenepy SS. D. Hedrick in October, 1967 to fit a Sinusoidal curve to LDV * 
* velocity data generated by an oscillating airfoil. The subroutine, utiliz- * 
* ing FOURIER, adapted from June 1977 OCS Airfoil program, fits a single -f 
* cycle sine wave to the data and computes the u- and v-velocity components . 
* at a phase angle specified for the analysis. > 
t ® 
i ® 


f®eeeeeneetretenee te teanhtetneaetret eee ee tenet eakeeeneteeee tee teeta ete eee eee ee tee e eee eeenenee 


IMPLICIT NONE 
CHARACTER*1 ANS, DUMMY 
CHARACTER*10 INFILE 
INTEGER NSAMPLES,I,PHASA,NSA,REJ,ARG 
REAL FU(0:3600),FV(0:3600),AU,AV,BU,BV,PI,X(3600),Y( 3600) 
REAL UVL,VVL,MENU,MENV, PREJ,MEANU(18),MEANV(18),MPHAVG(18) 
DATA PI/3.1415927/ 
WRITE(*,*)’DESIRE A VELOCITY/PHASE PLOT OF THE REFINED DATA?’ 
WRITE(*,*)’(Y/N--DEFeyY)’ 
READ(*,’(Al)}’) ANS 
TF (ANS 2EQ. °N’ .OR.-ANS 8EQ. ‘'n’) GO TO 7 
CALL RPLOT(NSAMPLES,INFILE,REJ, PREJ) 
7 OPEN(UNITel1l, FILE*INFILE,STATUS#’ OLD’ ) 
DO S I#l1,NSAMPLES 
READ( 117,100) X(1I),Y¥(T) 
S CONTINUE 
100 FORMAT(5X,F9.4,1X,F9.4) 
CLOSE(UNITe#11) 
DO 10 I=#1,NSAMPLES 
FU(I-1)=#X(I) 
FV(I-1)<«Y(I) 
10 CONTINUE 
NSA#=NSAMPLES-1 
CALL FOURIER(FU,NSA,AU,BU,1) 
CALL FOURIER(FV,NSA,AV,BV,1) 
DO 20 I#0,NSA 
FU(I+1)eMENU+AU*COS((2.*PI/FLOAT(NSAMPLES))*FLOAT(I))+ 


1 BU*SIN((2.*PI/FLOAT(NSAMPLES) )*FLOAT(TI) ) 
FV(I41)eMENV+AV*COS((2.*PI/FLOAT(NSAMPLES))*FLOAT(I))+ 
1 BV*SIN((2.*PI/FLOAT(NSAMPLES) )*FLOAT(1I1)) 


20 CONTINUE 
WRITE(*,*)’PLOTS WILL NOW BE PROVIDED TO COMPARE SINUSOID WITH’, 
i) ‘ VELOCITY DATA. PRESS <CR> TO CONTINUE.’ 
READ(*,’(Al)’) DUMMY 
CALL PLOT(NSAMPLES, INFILE, ?PU,FV,MENU,MENV,MEANU, MEANV, 
1 MPHAVG ) 
ARGeNSAMPLES * ( FLOAT( PHASA) /360. ) 
WRITE(*,*)’ ARG e@’,ARG 
UVLeFU( ARG) 
VVLeFV( ARG) 
RETURN 
END 
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SUBROUTINE PLOT(NSAMPLES, INFILE,FU,FV,MENU,MENV,MEANU, 
Sh MEANV, MPHAVG ) 


RAAARAHRAARARARARAARAARAARARARARARARARARKRARARAARARAARARARARARAARARARARAARARARARRARARARARARRARARE 


® ® 


* Written by S. D. Hedrick in October, 1987 to plot a comparison of velocity * 
* data and a sinusoidal curve fit. * 
* 
RARRAARARAARARAHRAARRARARARARARARARAARARRARHRARKAARKRRAAARRAARARARARARARARAARAARRARARARRARARRARE 
IMPLICIT NONE 
INTEGER NSAMPLES,1I,BOUND,STOPU,STOPV 
REAL *4 XVAL( 4000), YVALU( 4000) ,YVALV(4000) ,XCONTR(4),YCONTU( 4) 
REAL *4 YCONTV(4),FU(NSAMPLES), FV(NSAMPLES) 
REAL* 4 YMINU, YMAXU, YMINV, YMAXV, MENU, MENV 
REAL MEANU(18),MEANV(18) ,MPHAVG(18) 
CHARACTER*10 INFILE 
CHARACTER*11 XLABEL 
CHARACTER*15 YLABEL 
CHARACTER*21 UTITLE 
CHARACTER*19 VTITLE 
CHARACTER*1] DUMMY 
OPEN(UNIT#11,FILE=INFILE,STATUS#’OLD’ ) 
DO 10 Iel,NSAMPLES 
READ(11,100) YVALU(I),YVALV(I),XVAL(I) 
XVAL(I)@XVAL(1I)/10. 
10 CONTINUE 
100 FORMAT(5X,F9.4,1X,F9.4,8X,F4.0) 
XLABEL#=’PHASE ANGLE’ 
YLABEL@’VELOCITY (ft/s) ’ 
UTITLE#=’HORIZONTAL VELOCITIES’ 
VTITLE=’VERTICAL VELOCITIES’ 
YMINU=MENU-150. 
YMAXU@MENU+150. 
YMINVeMENV-15S0. 
YMAXVeMENV+150. 
XCONTR(1)=#6.0 
XCONTR(2)#0.0 
XCONTR( 3)=360.0 
XCONTR(4)=45.0 
YCONTU(1)=#6.0 
YCONTU( 2) =#YMINU 
YCONTU( 3) =YMAXU 
YCONTU( 4)=50.0 
YCONTV(1)=#6.0 
YCONTV(2)#YMINV 
YCONTV (3) =YMAXV 
YCONTV(4)=#50.0 
CALL LGPSPLOT(1,’EXSY’ ,XVAL,FU,NSAMPLES,XLABEL,YLABEL, , , , 
1 XCONTR,YCONTU, ,UTITLE) 
CALL LGPSPOINT(1,MPHAVG,MEANU,.03,3,18, , , ) 
WRITE(*,*)’PRESS <CR> TO TERMINATE PLOT’ 
READ(*,200) DUMMY 
CALL LGPS$TERMINATE PLOT(1,1) 
CALL LGPSPLOT(1,’EXSY’,XVAL,FPV,NSAMPLES,XLABEL,YLABEL, , , , 
XCONTR,YCONTV, ,VTITLE) 
CALL LGPSPOINT(1,MPHAVG,MEANV,.03,3,18, , , ) 
WRITE(*,*) ‘Press <CR> to terminate plot’ 
READ(*,200) DUMMY 
200 PORMAT(A1) 
CALL LGPSTERMINATE_PLOT(1,1) 
CLOSE(UNITe11) 
RETURN 
END 


td 


68 


SUBROUTINE PROPLOT(UNDIM,VNDIM,ZPO, XPOS,N) 
IMPLICIT NONE 

REAL*4 UNDIM(30),VNDIM(30),ZPO( 30) ,XCONTU( 4) ,XCONTV( 4) 
REAL*4 YCONTR( 4),UNDI(3000),VNDI(3000),2ZP(3000) 

REAL XPOS 

INTEGER N,NBR 

CHARACTER*29 XLABEU,XLABEV 

CHARACTER*32 YLABEL 

CHARACTER*27 UTITLE 

CHARACTER*25 VTITLE 

CHARACTER*1 DUMMY,ANS 

CHARACTER*9 XPO 

CHARACTER*4 XXPO 

XLABEUs’VELOCITY RATIO (U/U(INFIN))’ 
XLABEVe'VELOCITY RATIO (V/U(INFIN))’ 

YLABEL#’ PROBE POSITION (IN. FROM CENTER)’ 

UTITLEe’ HORIZONTAL VELOCITY PROFILE’ 

VTITLE@' VERTICAL VELOCITY PROFILE’ 

XPOS=XPO0S/3.0 

WRITE(UNIT#XXPO, FMTe#’(F4.1)'’)XPOS 

XPO#’X/C @'//XXPO 

XCONTU(1)#6.0 

XCONTU(2)20.0 

XCONTU(3)e1.2 

XCONTU(4)2#0.2 

XCONTV(1)2#6.0 

XCONTV(2)e-0.7 

XCONTV(3)=#0.7 

XCONTV(4)2#0.2 - 

YCONTR(1)=#6.0 
YCONTR(2)#-2.0 
YCONTR(3)=2.0 
YCONTR(4)#0.5 
NBR#10*(N-1)+1 
CALL LGPSPLOT(1,’EXSY‘',UNDIM,ZPO,0,XLABEU, YLABEL, 
i MP XGONTU, YCONTR, ,UTITLE) 

CALL LGPSPOINT(1,UNDIM,ZPO,.03,3.N, ,, ) 

CMLL LGPSPUT TEXT(1,.4,5.5,XPO, , ) 

WRITE(*,*)'PRESS <CR> TO TERMINATE PLOT’ 

READ(*,’(Al)’) DUMMY 

CALL LGPSTERMINATE PLOT(1,1) 

CALL LGPSPLOT(1,’EXSY’,VNDIM,ZPO,0,XLABEV,YLABEL, 

1 Poe XEONTV, YEONTR, ,VTITLE) 

CALL LGPSPOINT(1,VNDIM,ZPO,.03,3,N, ,, ) 

GRU LGPSPUTETEXT( 1 W4,5.5,XPO, , ) 

WRITE(*,*)’PRESS <CR> TO TERMINATE PLOT’ 

READ(*,’(Al)') DUMMY 

CALL LGPSTERMINATE PLOT(1,1) 

WRITE(*,*)’DO YOU DESIRE A SPLINE INTERPOLATION CURVE FIT OF’, 
1 ‘ THE DATA POINTS? (Y/N)’ 

READ(*,’(Al)’) ANS 

IF (ANS .EQ. ’N’ .OR. ANS .EQ. ‘n’) GO TO 10 

CALL LGPSPLOT(1,’EXSY’,UNDIM,ZPO,0,XLABEU, YLABEL, 

1 , » »XCONTU,YCONTR, ,UTITLE) 

CALL LGPSPOINT(1,UNDIM,2ZP0O,.03,3,N, +, ) 

CALL LGPSSPLINE QHC(UNDIM,ZPO,N,10,UNDI,2ZP, ) 

CALL LGPSPLOTM(1,’EXSY’,UNDI,ZP,NBR, ) 

CALL LGPSPUT TEXT(1,.4,5.5,XPO, , ) 

WRITE(*,*)’ PRESS <CR> TO TERMINATE PLOT’ 

READ(*,’(Al)’) DUMMY 

CALL LGPSTERMINATE PLOT(1,1) 

CALL LGPSPLOT(1,’EXSY’,VNDIM,ZPO,0,XLABEV, YLABEL, 
1 pu, XCONTV, YCONTR, ,VTITUE? 

CALL LGPSPOINT(1,VNDIM,ZPO,.03,3,N, ,, ) 

CALL LGPSSPLINE QHC(VNDIM,2PO,N,10,VNDI,2ZP,) 

CALL LGPSPLOTM(I,’ EXSY’ ,VNDI,ZP,NBR, ) 

CADEOEGPSPUT DEXT(1,.¢,5.5,XP0, , }) 

WRITE(*,*)’ PRESS <CR> TO TERMINATE PLOT’ 

READ(*,’(Al)’) DUMMY 

CALL LGPSTERMINATE PLOT(1,1) 

RETURN = 

END 


rf 
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